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Abstract 



We examine the reduced density matrix of the centre of mass on 
position basis considering a one-dimensional system of non- 
interacting distinguishable particles in a infinitely deep square 
potential well. We find a class of pure states of the system for 
which the off-diagonal elements of the matrix above go to zero as 
N increases. This property holds too for the state vectors which 
are factorized in the single particle wave functions. In this last 
case, if the average energy of each particle is less than a common 
bound, the diagonal elements are distributed according to the 
normal law with a mean square deviation which becomes smaller 
and smaller as increases towards infinity. Therefore when the 
state vectors are of the type considered we cannot experience 
spatial superpositions of the centre of mass and we may conclude 
that position is a preferred basis for the collective variable. 
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1. INTRODUCTION 



The old problem of how to derive the classical behaviour of macroscopic 
bodies from quantum principles now seems much closer to solution thanks 
to spontaneous localisation theories of the GRW type[|l| or by reference to 
decoherence which occurs for a system interacting with an environment [@] . 

In the case of a macroscopic body, decoherence is believed to suppress the 
off-diagonal elements of the spatial reduced matrix of the collective variables, 
such as the centre of mass, irrespective of initial conditions. This is equivalent 
to saying that we cannot experience spatial macrosuperpositions of the centre 
of mass and this explains why we do not encounter states of this kind when 
looking at everyday objects. Moreover, as long as we limit ourselves to 
considering spatial distributions of the collective variable and the observable 
quantities linked to this, we can interpret the reduced density matrix after 
dechoerence in terms of classical probabilities. 

Clearly if pure states exist for a closed macroscopic system, corresponding 
to a reduced density matrix for the centre of mass that is diagonal on the 
position basis, they will give the same conclusions for the description of the 
collective behaviour. This has led us to identify possible pure states of a 
system of many particles in an external potential with the above graceful 
property. We wish to make it clear that our purpose is this, not the study 
of a decoherence process. Decoherence touches us only insofar as the result 
it gives have provided the inspiration for this note. 

For our purposes we require a very simple model as described below. We 
believe that the final result will be to add something to a genuinely interesting 
subject, rather than just to present a textbook exercise. 

2. A VERY SIMPLE MODEL 

The model consists of a closed one-dimensional system of noninter- 
acting distinguishable particles with position coordinates xi, X2, ■ ■ ■ , xat in a 
infinitely deep square potential well V{x). This is actually a one-dimensional 
ideal gas in a container. The value of V{x) at the bottom of the well is taken 
as the origin of the energy scale. This region of zero potential occupies the 
interval (0, L) of the it is bounded on both sides by infinite potential 
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barriers. 

We will focus our attention on the above case, but will also consider 
infinitely deep potential wells U{x) represented by an analytical and posi- 
tive definite function which approach the form of V{x) in some limits. One 
candidate for U {x) could be: 

C/(x) = [2(x-L/2)/L]^- (1) 

with m a positive integer as large as we like it to be. 

Apart from the major simplification of no interaction between particles, 
a further simplification we adopt is to consider the particles as spinless and 
distinguishable. Each single particle must be identified by its own mass to 
achieve this, all the masses being different. However, this would introduce a 
slight but nevertheless unnecessary complication, and we therefore consider 
particles as having the same mass /j, but being of different colours: the visible 
spectrum does actually offer a large number of possibilities. 

Let us consider the energy eigenvalues problem for a single particle in the 
infinitely deep square potential well. The simple and very known result is 
the set of eigenfunctions: 



^nix) = \I'2/L smin-Kx/L) (n = 1, 2, ■ ■ ■ , oo) (2) 

which vanishes at extremes and L and are zero outside the well. These 
correspond to the discrete energy values: 

En = imrh)yi2fiL^) (3) 

and define a vector space for each individual component. A generic state of 
the system can be expressed using the tensor product space as: 

'^{xi,X2,- ■ • ,Xn) = Cm,n2,-,nNfm{Xl)(fn2{x2)---fnN{^N) ■ (4) 

Below we consider the factorized state vectors: 

N 

^{xi,X2,- ■ ■ ,xn) = '[[iJi{Xi) (5) 

i=l 

with reference to the system subject to the potential U{x) and require that 
the average energy of each particle is limited by a common bound W: 

2 2 

< iJi{xi)\ |- + U{xi) \^i{xi) >= ^|^+ < U{xi) ><W (6) 
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where pi is the momentum operator of the i particle. Since the average 
potential energy is positive, it follows that < > < 2iiW and clearly the 
mean square deviation of the momentum < pf > — < Pi >^< IjiW . Now 
indicating with (Tj =< xf > — < Xi the mean square deviation for the 
position of the i*^ particle the uncertainty relation imply that ai > / (S/iVT) 
for each i. This will be valid still in the limiting case of infinitely deep square 
well. 



3. THE REDUCED MATRIX FOR THE CENTRE OF MASS 

As we are interested in the collective behaviour of the system, we will use 
the position co-ordinate of the centre of mass: 

X = {xl + X2^ Vxn)/N 

and the coordinates of the particles with respect to their centre of mass: 
Xi — X -\- The are not independent, since they must satisfy the relation 
Hf=i 6 = 0. We will consider as independent variables X together with the 
first — 1 relative co-ordinates, i.e. 6, 6, ■ ■ ■ , ^Ar_i. The last relative co- 
ordinate is expressed as = — ii- The Jacobian determinant of the 
transformation is equal to N . We will write the integration volume element 
for the new variables as dV — dVc dVi, , where dVc — dX refers to the centre of 
mass and dVb — Nd^i d^2 ■ ■ ■ d^N-i to the internal degrees of freedom. The 
trace operations in which the continuous spatial reduced matrix of the centre 
of mass may be involved will therefore be executed using the integration 
element dX. 

We now consider a state vector of the system, ^(xi, ^2, • • • , xtv)- Its 
projection < X'|\E' > on a position eigenstate \X' > of the centre of mass is 
simply ^(X + ^i, X + ^2, ■ ■ ■ , ^ + ^n-i, X - J2f=i^ 6)- The reduced density 
matrix elements we are interested in are obtained by taking the partial trace 
of<X'|^>< ^\X > over the microscopic degrees of freedom: 

N-l N-1 

Qx'x^ / + E6)^(^' + 6'---'^'- E6)^H (7) 
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or taking the relation Y^iLi = into account with the aid of a delta function: 



N 

Qx'x = N 6(£ ii)^>*{XHu ■ ■ ■ , XHn)-^{X'Hi, ■ ■ ■ , X'Hn) d^i--- d^N 

i=l 



(8) 



4. THE OFF DIAGONAL ELEMENTS 



We will now look at the off diagonal elements gx'x- In the first place, we 
note that these are limited to the range (0, L) of X and X'. This is clear from 
observing Eq. and remembering that the wavefunctions of the system are 
zero when the position coordinates Xi, X2, ■ ■ ■ ,X]y are outside this range. We 
may also note that, for the same reason, only the integration region for which 
both inequalities —X < < L — X and —X' < < L — X' are satisfied 
for each i may contribute to the integral in Eq. which defines the matrix 
elements. Let X' > X and A = X' — X. Therefore, indicating as f2 the 
volume spanned by the values of the in the intervals —X < < L — X' 
and using the expression in Eq. for the generic state vector: 

QX'X = X C„i,---,njv ^rni,---,mM -^ni,---,njv;»ni---,mjv (9) 

where 

N N N 

In.,-,n^;m„-,m^ = J Si^. ^ ^^X + ^ ^^^{X' + . (10) 

i=l 1=1 j=l 

Applying the Schwartz inequality: 

N N N 

i=l 1=1 j=l 

(11) 

Since the term to be integrated in the first factor is a definite positive quan- 
tity, we can remove the delta function which restricts the integration region 
and obtain a fortiori: 



(12) 



i=l " j=l ■ 
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If the integration variables are changed the previous expression may be 
rewritten as: 



rL-A N L 

i=l-'^ 3=1 



Now if we use 5a to indicate the maximum value with respect to the index 
n of the quantity: 

j-L-A i-L i-A 

/ \ipn{x)\^ dx = I \ipn{x)\^ dx = I - j \ipri{x)\^ dx 

Jo J A Jo 

where we used the simmetry properties of the ipn{x) and the fact that the 
same wavefunctions are normalized, we finally have: 

|-^n.i,---,njv;'ni,---,r?ijv I — "^A • (■^'^) 

To return to the Eq. (|^) above and using the triangular inequality we finally 
obtain: 

\0xx'\<NS5^ (15) 
where S indicates the double sum: 

|Cni,---,njv I |c?ni,---,mjv I — ( E/ l'^rii,---,njv I) 

Quantity 6a is clearly less than one for each finite value of A. Therefore, 
if the value of S is given by a polynomial in N, the non-diagonal matrix 
elements tend to zero in the limit of large A^. The results obtained refer 
to those matrix elements for which X' > X, but may also be extended 
immediately to the other half, remembering that our matrix is hermitean. 

If only a number / of terms appear in the expression of the generic state 
vector, one can easily see that, given the normalisation constraint, the max- 
imum possible value for S is /. Therefore, the condition required to have a 
corresponding spatial reduced density matrix for the centre of mass which is 
diagonal is clearly satisfied by the fundamental state, by the low-lying energy 
eigenstates and by their superpositions. 

In the case of a factorized state vector as in Eq. (^, by applying the 
Schwartz inequality once again and after some calculations like above, we 
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obtain: 

\qx'x\' <N^f[(l- r mx,)\'dx)j f[(l- ['^\^,{x,)\'dx, 

i=l \ JL-A ) -^-y \ Jo 

(16) 

If we exclude very unusual forms of the functions, all factors on the 

right-hand side are less than unity and their product therefore decreases 
rapidly as N increases. To be more precise, we assume that the particles 
have the same individual wave function. In which case: 

\g^,^\<Narbr (17) 

with ttA, &A positive quantities less than one for any finite A and the off- 
diagonal matrix elements exponentially tending to zero as N increases. 



5. THE DIAGONAL ELEMENTS FOR FACTORIZED STATES 

We will now examine the diagonal elements of the reduced matrix, with 
the reference only to a factorized state vector as indicated in Eq. (|), as this 
is the only case in which we can obtain specific results. Using Eq. (H) and 
the Fourier representation of the delta function, we have: 

N r ^ /■ 

Qxx = — dkH e^'^HnX + ^,)ij,{X + e.) d^, (18) 

which may be rewritten as: 

Qxx = ^j rfA;e-*'=^^ \{ j e^^^'^;(a;,)^.(x.) dx, . (19) 

Apart from the factor N (since we are using the differential element dX 
instead of dY here), this is nothing other than the expression for the prob- 
ability density distribution of the sum NX = Y = xi-|-a:2-l-... + a;Ar. We 
could have written this directly using the characteristic functions method by 
observing that the Xi are distributed in accordance with the independent 
probability densities Fi{xi) = ip*{xi)ipi{xi). 

We will now consider the second central moments cxj =< (xj— < Xj >)^ > 
and the sum sn = YJi=i If we require that the average energy < Ei ><W 
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for each particle, it follows as noted that the second central moments are 
limited by a common bound independent of N. This means cxj > a were a 
is a positive constant. 

In the case of the infinitely deep square well, since the F, carried 
by a finite interval, the only requirement that each cTj > a is sufficient to 
satisfy the Lindeberg conditionP] for the Central Limit Theorem (CLT). 
Therefore, under the conditions indicated, the distribution of the diagonal 
elements gxx will be well expressed for a very large by the Gaussian: 



where < a > indicates the quantity sn/N and < X >= < Xi > /N 

is the average value of the centre of mass position. < X > and < a > will 
depend generally on time, but the Gaussian form for large will in any case 
be guaranteed since the conditions for the CLT are satisfied independently of 
time. Moreover, as the wave functions are zero outside the interval (0, L), the 
particles mean square deviations are also limited from above. We certainly 
have CTj < for each particle. Therefore < a >= J2f=i <^i/N remains limited 
and < a > /N goes to zero as increase. On the time dependence of the 
central value < X > and of the variance it can be said that this dependence 
will be of the periodic type, since the Bohr frequencies involved are whole 
multiples of a basic frequency, as can be seen from Eq. (^) 

In the case of the potential U{x), we once again find that the distribution 
of diagonal matrix elements follows the normal law if in addition to the 
requirement that the average energy of each particle is limited (which means 
that their second central moments are limited as well) we also consider the 
third absolute central moments z/j =< | Xj— < Xi > \^ > and make the 
physically reasonable assumption]^ that Vi < (3 for some positive constant 
j3. Indeed in this case the Liapunov conditions for the CLM|^ are satisfied. 



We have considered the reduced density matrix of the centre of mass on 
position basis corresponding to the pure states of a one-dimensional system 
of non-interacting distinguishable particles in an infinitely deep square well 




(20) 



6. FINAL REMARKS 
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potential. By expressing the generic state vector in terms of energy eigen- 
states, a sufficient condition was identified for the matrix in question to be 
diagonal The condition is that the quantity S defined in the text is hmited 
by a polynomial in the number of particles. This is clearly satisfied by 
the fundamental state, by the low-lying energy eigenstates and by their su- 
perpositions. Therefore, when the system is in one of these pure states, or 
others for which the condition on S is verified, we cannot experience spatial 
superpositions of the centre of mass, its spatial reduced density matrix can 
be thought of in terms of classical probability and to conclude the position 
will appear as a preferred basis for the collective variable. However we are 
cautious about this last conclusion, since non-diagonal elements of "small" 
subsystems are known to be small in an arbitrarily chosen basis. There even 
exist theorems for the most probable value of entropy of such systems]^]. 

Since we expressed the generic state vector by means of energy eigenstates 
the coefficients of the expansion depend on time only by way of a phase fac- 
tor, and as S depends on their moduli when the condition required for this 
quantity is satisfied, it is satisfied permanently. The above states for our 
system of non-interacting particles do not therefore appear to be the final 
result of a decoherence process. Nevertheless, if we introduce an interaction 
between the particles ( implying a time dependence for the moduli of our 
expansion coefficients ) we could think of the possible states which evolve 
fulffiling at the end the condition on S as the result of a process of deco- 
herence, considering our gas as consisting formally of a "collective system" 
described by the centre of mass and an ( internal ) environment described 
by the microscopic degrees of freedom. Coupling between the two formal 
systems is guaranteed by the external potential. 

We also considered the case in which the system's state vector is factor- 
ized in the individual wave functions of its components. This state vector 
is a long way from being generic, but this assumption is natural if there is 
no interparticle interaction. Here too there are indications that the corre- 
sponding spatial reduced density matrix of the centre of mass is diagonal. 
Moreover, if we require the average energy of each particle to be less than 
a common value, the conditions for the CLT are satisfied and the diagonal 
elements of the matrix considered are well expressed by a Gaussian distri- 
bution as increase. The Gaussian parameters are time dependent, but 
the mean square deviation decreases in any case as This is consistent 

with a recent work in which we discussed the localisation properties of a 
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macroscopic body with low total average energy |^. Factorised states of such 
a type appear for any practical purpose as localised states of the centre of 
mass. Each of them appears potentially to be the result of successive spon- 
taneous random localisations of the individual particles, as contemplated by 
theories of the GRW type. An appropriate classical probability distribution 
of similar states localised at various values of the position can at a given 
moment mimic the distribution of diagonal elements of the reduced matrix 
corresponding to the states considered before. But the meaning of the two 
cases is very different [|] . 

To conclude, we must remember that some time ago, the example of a 
macroscopic ball bouncing elastically between two parallel walls positioned 
normally to the x axis was used to demonstrate the difficulties encountered 
in describing a similar situation with quantum mechanics[§]. The gas in our 
example cannot be considered a compact body like a ball, but we agree that 
it is macroscopic. If this is described by a factorised state vector like those 
indicated the centre of mass position is distributed according to a well lo- 
calised Gaussian, whose central value moves with periodic motion between 
the two walls. However we do not claim that this distribution, which is a 
general consequence of the statistical independence of the particles, repre- 
sents classical behaviour which would require a Gaussian distribution about 
a time-dependent centre for each wave function in the ensemble. This could 
be expected only for a bound state describing a macroscopic object, not for 
a gas. 
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